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In magnetic fields applied within the angular range of the surface superconductivity state a mag- 
netically anisotropic layered medium is created in structurally isotropic, sufficiently thick niobium 
films. Surface (Kulik) vortices residing in the superconducting sheaths on both main film surfaces in 
tilted fields are shown to undergo a decoupling transition from a coherent to an independent behav- 
ior, similar to the behavior observed for Giaever transformer. At the transition a feature in pinning 
properties is measured, which implies different pinning for the lattice of surface vortices coherently 
coupled through the normal layer and for two decoupled vortex arrays in the superconducting surface 
sheaths. 

PACS numbers: 74.25.Op, 74.25.Qt, 74.25.Ha, 74.78Db 



I. INTRODUCTION 

The existence of vortices and vortex states within the 
surface superconductivity (SC) state has been discussed 
since the discovery! of this phenomenon within the ap- 
plied magnetic field range B c2 < B a < B c3 (B c2 and 
B C 3 are the second and surface superconductivity critical 
fields). The existence of Abrikosov vortices^, Abrikosov- 
like stated, giant (multi-quantum) vortex stat o 5 i 6 i 7 , and 
Kulik (surface) vortices&^iiS above B c2 were suggested 
depending on sample size and shape and applied field 
orientation. In this work, we mainly deal with the mag- 
netic behavior governed by surface vortices in the su- 
perconducting surface sheath. The structure of a sur- 
face vortex, to a large extent, reproduces the structure 
of an Abrikosov vortex with the length equal to the 
thickness of the superconducting layer (surface sheath). 
Fink and Kessinger showecUi that the thickness of this 
sheath d sc ~ 1.6£(T) at B a ~ B c2 and approaches £ at 
B a = B C 3 for a superconductor with k = 10 (re = A/£ 
is the Ginzburg-Landau parameter, £ and A are the co- 
herence length and magnetic field penetration depth) . In 
the case of a relatively thin superconductor, as the films 
investigated in this work, the superconducting surface 
sheaths on the main film surfaces are separated by a nor- 
mal layer 

d n ~d p - 3.2£(T) (1) 

with dp being the film thickness. 

For a certain magnetic field range applied at an an- 
gle (9) to the surfaces, two independent flux-line lat- 
tices (FLLs) can be formed in thin films and layered 
systems ASiSiiSiiSiiiiiiSiiSiii This co-existence is possi- 
ble due to strong structural or magnetic anisotropy and 
to the two components of the applied field. §*i^ In this 
case, the out-of-plane field component (B a ±) would be 
responsible for the out-of-plane (perpendicular) FLL and 
the in-plane (parallel) component (-B a ||) for the in-plane 



vortex lattice. The co-existence of two FLLs has experi- 
mentally been shown for structurally isotropic films simi- 
lar to those investigated in this work^i^ii^i&. Therefore, 
we will hereafter assume that in the films investigated two 
co-existing FLLs, perpendicular and parallel to the main 
film surface, are present at fields tilted with respect to 
the surface. Fig. da) schematically shows the two-FLLs 
structure. 

We further assume that at B a > B c2 and over the an- 
gular range of \9\ < 40°, which is the range over which 
the surface SC state can be measured*^ (i) the paral- 
lel FLL in thick films d p » t; transforms into a giant 
vortex above £? c2 £Li£ii£ This is a reasonable assump- 
tion because the surface SC state, just like the in-plane 
Abrikosov vortex rows, forms due to the -B a | compo- 
nent. The shielding super-currents of the giant vortex 
flow within the surface sheath, (ii) The perpendicular 
FLL forms two arrays of quasi-2D surface vortices resid- 
ing in the superconducting sheaths on both main sur- 
faces of the films' (Fig. ^b)). The co-existence of the 
giant vortex and the surface vortices was discussed in 
Refs. IglfToL It was also shown that Kulik vortices can 
form triangular or square lattice, depending on the ap- 
plied field and its orientation 9& 

Neither theoretical nor experimental work has been 
published, which shows any kind of interaction between 
the two surface vortex arrays on the opposite surfaces of 
a thin flat sample in the surface SC state. We expect that 
the arrays can behave either coherently or independently, 
depending on d n , d sc , B a and 9. These parameters can 
affect the coupling force between the vortices as it was 
shown for the case of two magnetically coupled supercon- 
ducting films (superconducting Giaever transformer—) in 
fields B a < B c2 m2£& At the transition between the co- 
herent and independent regimes a small pinning change 
would be expected. Therefore, one needs a technique 
sensitive to such small changes in fields applied nearly 
parallel to the film surface. Such fields are necessary to 
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FIG. 1: The vortex arrangements are schematically shown for 
the film cross-section in a field B a applied at an angle 9 to the 
main film surfaces (a) below B C 2, and (b) above B C 2 after the 
decoupling transition. The black layer near the circumference 
of the film denotes the surface SC sheath with the light gray 
surface vortices in it. The black circles in (a) imply in-plane 
Abrikosov vortex rows-*- parallel to the surfaces. The dark 
grey stripes perpendicular to the surfaces show the out-of- 
plane Abrikosov vortices. Although it is not clear on the 
picture but it is assumed that vortex lines of both Abrikosov 
lattices in (a) do not cross one another. 



enable the surface SC state. In this work, we describe re- 
sults of mechano-magnetic experiments on niobium (Nb) 
films of different thicknesses and provide experimental 
evidence for the decoupling transition of the surface vor- 
tices within the surface SC state. 



II. EXPERIMENTAL DETAILS 

The experimental technique employed in this work is 
based on mechanical oscillations of a superconductor at- 
tached to a Vibrating Reed (VR) in an external magnetic 
field (for a review see Ref. 21 and references therein). 
The experimental VR setup employed in this work can 
be found in Ref. 0. This technique is very sensitive to 
magnetic properties of the superconductor, in particu- 
lar to the pinning of vortices, in fields parallel to its 
largest surface (9 = 0°). The physical reason for this 
sensitivity is that the very small field component ~ B a ip 
perpendicular to the applied field, arising when the su- 
perconductor is tilted by a very small angle (typically 
(p < ICT 5 degrees), is shielded by superconducting cur- 
rents (generally defined by the pinning of vortices). The 
shielding currents cause the external field to curve around 
the tilted superconductor. This field distortion leads to 
an additional line tension (stiffness) in the system, which 



is proportional to the increase of the length of the field 
lines near the edge of the tilted superconductor. As a 
result, the resonance frequency (ui) of the VR with the 
attached superconducting sample increases with field. If 
the shielding supercurrents become smaller, for example 
in the vicinity of the upper critical field or the critical 
temperature (in general due to a pinning reduction), the 
resonance frequency decreases. The damping (r) of the 
oscillator, which is measured simultaneously with the res- 
onance frequency, is proportional to the corresponding 
energy dissipation occurring in the oscillator (reed plus 
superconducting sample) due to vortex movement and 
the internal friction of the reed material. The peak in 
the damping, usually measured as a function of field or 
temperature, corresponds to the vortex depinning line. 22 
However, in the case of the surface superconductivity (the 
so-called giant vortex state) the peak can have a differ- 
ent origin related to the shielding/pinning properties of 
the giant vortex (see, for example, Ref. |y and references 
therein) . 

In the case of a thin conventional superconductor, such 
as Nb-film, the behavior of the VR as a function of tem- 
perature and field is, to a large extent, governed by the 
magnetic properties of the surface of the superconduc- 
tor. Therefore, the unique properties of the VR tech- 
nique should allow us to detect changes in the behavior 
of the surface vortex arrays, whose pinning can influence 
the shielding property. Accordingly, the resonance fre- 
quency change oj 2 (B a ) — w 2 (0) and the damping T(B a ) 
measured in the experiment are respectively expected to 
provide information on surface vortex pinning and energy 
dissipation produced by vortex movement. 

The increase in the resonance frequency vanishes as 
soon as vortex pinning and the shielding become negli- 
gible. Thus, we can measure not only B C 2 at 6 — * 90° 
and B c3 at 9 — ► 0°, but also the angular dependence of 
the upper critical field (j? MC (fl))' 3 ' 4 ' 15 Naturally, we define 
B uc (9 = 0°) = B c3 and B uc (9 = 90°) = B c2 . 

Another important feature of our experimental setup 
is the high angular resolution of the rotation system^ As 
a consequence, the angle 9 of the field with respect to 
the main film surface was defined with an accuracy of 
±0.01° in the vicinity of 0° and of ±0.5° at 6 > 3°. At 
9 S> 0°, the angular resolution is smaller for these exper- 
iments due to a small A 2 B CU /(A9) 2 at 9 ^ 0° (see Fig. 5 
in Ref. 0). The accuracy near 0° is limited by the sensi- 
tivity of the Si-oscillator onto which the superconducting 
sample is attached. This Si-oscillator has a quality factor 
Q ~ 10 6 at the temperatures of the measurements. 

Nb-films of different thicknesses d p ~ 120 nm (Nbl20), 
400 nm (Nb400), and f200 nm (Nbf200) were in- 
vestigated in this work. The polycrystalline films 
were sputtered onto an oxidized silicon wafer at room 
temperature^ Superconducting properties of these films 
were characterized in earlier works4»^i& The coherence 
length at zero temperature £(0) for all the measured sam- 
ples was estimated to be ~ f2.5 nm (£(T = 5K) ~ 
19 nm). We estimate a Ginzburg-Landau parameter 
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FIG. 2: (a) The resonance frequency change Auj 2 = u 2 (B a ) — 
u> 2 (0), (b) the corresponding damping, (c) the first derivative 
of Auj 2 on B a with the arrows denoting the minima (B m ), and 
(d) the enlargement of the damping onset at Bonset which is 
marked by the arrows for the Nb400-film. The dashed and 
dashed-dotted lines mark B C 2 and B c s, respectively. 
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FIG. 3: The same as in Fig.EJfor the Nbl200-nlm. 
A/£ ~ 10 for the three films. 4 



III. EXPERIMENTAL RESULTS 

Figures |21 and |3| show the measured resonance fre- 
quency change oj 2 (B a ) — w 2 (0) (a) and the VR damp- 
ing T(B a ) (b) as a function of applied field at different 
angles 9 and at fixed temperatures for the Nb400 and 
Nbl200-films, respectively. The key-feature in these fig- 
ures is the appearance of an unusual non-monotonic be- 
havior at angles 9 > 9™, which is best seen as a minimum 
at a field we define as B m in the first derivative of the 
resonance frequency (Figs. |2Ic) andEJc)). This "criti- 
cal" angle 9™ is ~ 0.45° and 0.16° for the Nb400- and 
Nbl200-films, respectively. B m (9) for both films is shown 



in Fig. Ufa). The angular range of the non-monotonic 
behavior (9™ < 9 m < 20°) is within the range of the 
surface SC existence4^ In principle, one would tend to 
observe the position of the minima (B m ) at B C 2- In this 
case, the minima would naturally indicate a change in the 
shielding property when the bulk superconductivity col- 
lapses and only surface superconductivity persists. How- 
ever, the minima do not coincide with the experimentally 
measured B c2 (marked by the dashed lines in Figs. [21 El 
and^}. Instead, the i? m -behavior is more complex being 
angular dependent. We argue below that this behavior 
can be explained as decoupling of the coherent arrays of 
the surface vortices, which undergo a decoupling transi- 
tion from coherent behavior at B a < B m to independent 
behavior at B a > B m . This transition is promoted by 
the magnetically anisotropic medium created in the films 
in fields within B C 2 < B a < B c s applied nearly parallel 
to the surface. Taking into account the thickness of the 
films, the coupling between the coherent surface vortex 
pairs is of magnetic nature. 

IV. DISCUSSION 

A. Angle dependence of B m 

At 9 < 9™, B a ± is too small to induce a large density 
of surface vortices in the surface sheaths. This is likely to 
imply that the magnetic behavior in this angular range 
is overwhelmingly governed by the giant vortex surface 
shielding current. Therefore, the VR signal is not sensi- 
tive enough to reveal a possible decoupling transition for 
only a few surface vortices up to B cu . As 9 becomes larger 
(Fig. Ufa)), B a ±_ increases and more coupled surface vor- 
tices are created. These vortices start interacting within 
each sheath at a characteristic crossover field given bySP. 

Bcr = vktfr ' (2) 

where $o is the flux quantum and ag r is the intervortex 
distance at the crossover. As soon as intervortex interac- 
tion between surface vortices within one sheath becomes 
stronger than the coupling force, the decoupling between 
the coherent pairs of surface vortices takes place, form- 
ing two independent arrays of (2D-like) vortex lattices, 
one in each surface sheath. The magnetic coupling of the 
surface vortices is weak due to the relatively large dis- 
tance d n . Thus, relatively weak in-plane vortex- vortex 
interaction should be enough to decouple the coherent 
behavior. A sufficiently strong intervortex interaction 
for a decoupling would occur at an intervortex distance 
ag r ~ 2A(T). If we assume that the intervortex spacing 
for the triangular surface vortex lattice ism 

a Q ~ (2$ /\/35 a sin 9) - 5 , (3) 

one finds that at 9™, a Q ~ 0.5 /xm for the Nb400-film and 
oq ~ 1 fim for the Nbl200-film. Thus, the intervortex 
spacing is of the order of 2A(T = 5K) ~ 0.38 fim at 9™. 
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FIG. 4: (a) Normalized field B m at which minima are ob- 
served in the resonance frequency change, and (b) damping 
onset as a function of angle. The arrows mark 9™ in (a) co- 
inciding with the step-like feature in (b) for the Nb400- and 
Nbl200-films. The dash-dotted, nearly vertical line in (a) 
shows Bcr(0) obtained from Eq. assuming a fixed value 
for the intervortex distance ao = ag r = 2A(5 K) = 380 nm. 
The dotted line in (b) shows the plateau at small angles and 
an additional experimental point measured at 9 — 0° for the 
films, which cannot be shown in the logarithmic scale. 



Let us assume that the decoupling transition occurs 
when the perpendicular component of the applied field 
B a ± = B a sm9 is equal to B cr . Then, the dash-dotted 
line in Fig. E{a) shows B cr (8) oc (o^sinf?) -1 expected 
from Eq. © with the fixed a = a c r = 2A(T = 5K) = 
380 nm for the NM200- and Nb400-films. As can be 
seen, Eq. with the fixed ao does not describe the 
B m (9) behavior over the entire decoupling line except at 
the angle ~ 1.5° at which the calculated curve crosses 
the corresponding experimental i? m -lines. The disagree- 
ment between the experimental curves and Eq. J3J should 
actually be expected. Indeed, in the case of Eq. © 
only one parameter - ao(B a ±), - affecting the decou- 
pling, changes with field (or angle). Whereas in our case 
there are at least four variables affecting the coupling: 
ao(B a ±), dsc(B a \\), d n (d sc ) and surface vortex pinning. 
The behavior of the B m (9) curves in Fig.^Ja) can be ex- 
plained by four major factors, which influence the above 
mentioned variables: (i) the nature of the decoupling 



which occurs within the surface SC state and with an 
enormously large interlayer spacing d n and small d sc , (ii) 
surface roughness of the measured films, (iii) surface vor- 
tex pinning which is not accounted for in Eq. @, and 
(iv) angular dependence of the upper critical field B cu . 

The first factor responsible for the B m behavior below 
the crossing point is particularly well described on the 
example of the thicker film (Nbl200) with much weaker 
coupling (due to the larger d n ) than that in the Nb400- 
film. For the Nbl200-film the ^-independent plateau is 
observed for B m at 9 < 2°. In this range the decou- 
pling is driven by the reduction of d sc above B c i with 
increasing B a mL The smaller d sc (the larger d n ) leads to 
a reduction of the coupling and pinning for both surface 
vortices of all coupled pairs. As soon as the decoupling- 
threshold is reached, the surface vortices from each pair 
are likely to be dragged apart by shielding currents in- 
coherently oscillating on opposite surface sheaths. As 9 
approaches the crossing point at ~ B cr , B m starts to 
curve downwards being also affected by intervortex in- 
teractions. A similar, but stronger effect experiences the 
thinner film (Fig. Ufa)). The stronger angular depen- 
dence below the crossing point is likely observed because 
the decoupling threshold is higher than for the thicker 
film. Hence, to reach the decoupling higher fields have 
to be applied (Fig. Ufa)), which result in larger surface 
vortex populations and, consequently, in stronger inter- 
vortex interactions. 

The second factor can have some influence in the vicin- 
ity of the crossing point. The values of ao, calculated for 
"ideal" film surfaces, are likely to be underestimated due 
to the surface roughness present in real films ^ The sur- 
face roughness model implies that even if 9 = 0°, the 
flux would intercept some localized areas of the rough 
surface. Thus, the rougher the surface, the more surface 
vortices are expected to populate the sheaths in applied 
fields nearly parallel to the surface. The Nbl200-film was 
found to have a larger value of the root mean square sur- 
face roughness (6.6 nm) than the Nb400-film (5.3 nm)^ 
This result can contribute to the fact that the decou- 
pling has been observed starting from a smaller 9™ for 
the Nbl200-film than for the thinner film. Apparently, 
at larger angles the surface roughness factor becomes less 
significant. 

The third factor is responsible for the disagreement at 
larger angles. The pinning experienced by the surface 
vorticesi2iS4i2^, which we neglected to a large extent in 
the above consideration, is likely to modify significantly 
Eq. derived by assuming pinning-free environment, 8 

The fourth factor can also influence B rn in particular 
at larger angles, since, for example, B cu at 5° is about 
10% smaller than B c3 M^ This can affect the thickness of 
the surface sheath d sc , pinning and shielding properties, 
and, therefore, the decoupling. 

In Fig. 01a) one sees that increasing 9, B m approaches 
B C 2- It may seem surprising that the feature attributed 
to the decoupling in the surface SC state still exists 
below B C 2- However, it was shown in a number of 
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theoreticalSi^SiSi and experimental worksi^iLSS that a gi- 
ant vortex state within a superconducting surface sheath 
can be nucleated at sufficiently high fields below B C 2- In 
this case, the magnetic anisotropy (the layered structure) 
of the surface SC state can also be preserved below B C 2- 
In addition, the superconducting order parameter within 
the d„-layer is substantially reduced due to a large num- 
ber of densely packed in-plane Abrikosov vortices^i In 
this case the magnetic anisotropy (layered structure) is 
effectively maintained below B C 2- We stress that the de- 
coupling is observed only within the angular range of the 
surface SC state existence (0 < 9 < 40°)^i^ defined as 
B cu > B C 2- Therefore, the decoupling appears to be a 
realistic scenario below B C 2, as well. 

Summarizing, the decoupling behavior in Fig. 0] can 
be described as follows. At 9 < 9°[ l , the giant vortex 
shielding overwhelmingly dominates, so that the possi- 
ble decoupling of only very few surface vortices cannot 
be detected by the VR technique: neither the measured 
resonance frequency change nor damping show an un- 
usual behavior up to the vicinity of B C 3 as if there were 
no transition at B C 2- At 9 > 9™ the decoupling occurs 
at B m {9). At fields below the decoupling pinning and 
shielding properties behave in a usual way as described 
for the VR technique As the field further increases 
the coupling force between the sufficiently large amount 
of surface vortex pairs becomes too small to prevent the 
decoupling. The decoupling is likely driven by two dif- 
ferent mechanisms below and above the crossing point: 

(i) Below the crossing point, as the coupling forces be- 
come too small due to the d sc reduction with increasing 
field, the independent surface vortices become more mo- 
bile due to incoherent oscillation of the shielding currents 
on the opposite film's surfaces. As the result, the shield- 
ing properties slightly weakens (resonance frequency), 
and the dissipation notably onsets (damping). However, 
the mobility of the vortices is expected to be incomplete 
due to the arising pinning of individual surface vortices. 

(ii) Above the crossing point, this scenario is further com- 
plicated by an additional parameter: intervortex inter- 
action, which assists in the decoupling process. In this 
case, the mobility of the vortices would be restricted by 
a collective process which arises from interplay between 
vortex- vortex interaction and pinning of the surface vor- 
tices. 

In both cases, the shielding would be slightly weakened 
at the decoupling and partially regained after pinning in- 
dependent 2D-like surface vortex arrays in the supercon- 
ducting surface sheaths. Experimentally, this behavior 
has produced the observed minima in the first derivative 
of the resonant frequency change (Figs. Etc) and^Jc)) 
and the apparent enhancement of the damping at -Bonsot 
(see the arrows in Figs. |2^d) and[«Hd)). Note that in 
Fig-EIb), the B onset (9) dependences for the Nb400 and 
Nbl200-films show step-like features at 9°[ l . These steps 
coincide precisely with the appearance of B m in the res- 
onance frequency change for both films. 
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FIG. 5: The maximum coupling force (F C m) between two sur- 
face vortex arrays on the opposite film surfaces as a function 
of the intervortex distance (ao(B a ±)) between the surface vor- 
tices (Eq.^J. The dotted line indicates ao = 2A = 0.38 /im at 
T= 5 K. 



B. Coupling force 

For simplicity, assuming that the thickness of the 
superconducting surface layer d sc ~ 1.6£ (in fact, 
d sc (B a ±)—), the maximum coupling force (F cm ) at high- 

20 



flux density regime (B > B cr ) is given by 2 * 
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1 - exp( d sc ) 



exp( d r 

ao 



), 



(4) 

provided that d sc /\ -C 1, d n -C A 2 /d sc , and £ <C 
ao/(27r) <C A, where /io is the permeability of free space. 
However, we should note that some of the actual condi- 
tions for the Nbl200-film at T = 5 K are slightly softer 
than the given above: d n < X 2 /d sc and ao/(27r) < A. 
We believe that it should be acceptable for our estimate, 
especially taking into account that the main condition 
of thin superconducting layers d sc /X < 1 is fulfilled. 20 
In addition, we also note that A becomes larger and d sc 
smaller with increasing field, reinforcing the applicability 
of the corresponding inequalities. 

In Fig. [21 F cm as a function of ao(B a ±) is shown for 
all the films. At the decoupling crossover (9 m l ) indicated 
by the arrows, F cm ~ 8.9 x 1CT 16 N and 5.9 x 1CT 17 N 
for the Nb400- and Nbl200-films, respectively. 

Taking into account the trend for thinner films to pro- 
duce the decoupling onset (6™) at larger angles and 
stronger magnetic fields (Fig.^Ja)), the expected decou- 
pling onset for the Nbl20-film would be at or slightly 
below ao — 2A (the dotted line in Fig. |SJ). In this re- 
gion, F cm is much larger for the Nbl20-film than that 
for the thicker films. Importantly, it is nearly within the 
region where F cm is nearly independent on B. There- 
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fore, to observe the decoupling transition a much larger 
B a ± (smaller ao), implying larger 6, should be applied in 
order to reduce F cm and to increase the intervortex inter- 
action. However, a minimum in the resonance frequency 
corresponding to the decoupling was not observed for this 
film, nor the step-like feature in the behavior of B onse t(9) 
(Fig-HP 5 ))- Instead, B onsct (8) / B c2 is clearly larger than 
for the thicker films. This behavior indicates stronger 
shielding and pinning, which remain unaffected by the 
decoupling but affected by the critical field dependence 
B cu {6) only. 

V. CONCLUSION 

In summary, from VR experiments on thin Nb films 
we have obtained evidence indicating that for sufficiently 
thick films the surface vortices in the surface SC state 
undergo a decoupling transition. At small fields/angles 
the aligned vortices are coupled through the normal layer 
d n exhibiting a coherent 3D-like vortex lattice behav- 



ior. At larger fields/angles the surface vortices decouple 
forming two independent vortex arrays (2D-like behav- 
ior) in the superconducting surface sheaths. In films with 
dp ^> the coupling between the aligned surface vor- 
tices appears to be too strong so that the experimental 
observation is not possible with the VR technique. By 
comparison, we note that the loss of the 3D coherence 
in a lattice of aligned pancake vortices in layered high- 
temperature superconductors was explained in terms of a 
melting phase transition from the 3D vortex pinning state 
to the regime of independently pinned 2D vortex lattices 
for Josephsoni£ and magnetic couplings between the 
superconducting layers. 
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